Main Synchronization Sequence Design Method for Global Covering Multi-Beam Satellite LTE

ABSTRACT

Disclosed is a main synchronization sequence design method for global covering multi-beam satellite LTE, comprising the follow steps: extending a standard Zadoff-Chu sequence to a generalized Zadoff-Chu sequence so as to obtain an initial candidate main synchronization sequence set of more candidate sequences; gradually narrowing down the candidate main synchronization sequence set according to a selection standard of a main synchronization sequence to obtain a final candidate main synchronization sequence set; and obtaining a main synchronization sequence set with optimal eclectic performance and complexity from the final candidate main synchronization sequence set. According to the present invention, a main synchronization sequence with optimal eclectic performance and calculation complexity can be designed for a global covering same-frequency networking multi-beam satellite LTE system.

FIELD OF THE INVENTION

The present invention belongs to wideband wireless communication field, and particularly relates to a main synchronization sequence design method for global covering multi-beam satellite LTE.

BACKGROUND OF THE INVENTION

Satellite mobile communication is one of the necessary means for communication anywhere. In recent years, the fourth generation (4G) terrestrial cellular mobile communication became matured gradually, and the 3rd Generation Partnership Project Long Term Evolution (3GPP-LTE), which is based on key techniques such as Orthogonal Frequency Division Multiplexing (OFDM), Multiple-Input Multiple-Output (MIMO), and same-frequency networking, etc., has been commercially available gradually. How to apply Terrestrial LTE (T-LTE) that has characteristics such as high speed, high capacity, high spectrum efficiency, and high power efficiency, etc. in satellite mobile communication to establish multi-beam Satellite LTE (S-LTE) mobile communication systems based on same-frequency networking technique is a hotspot and a challenge in the current research in satellite mobile communication field.

Providing mobile communication services with Geostationary Earth Orbit (GEO) satellites has many advantages: theoretically, a single GEO satellite can cover 42.2% earth surface, and 3 GEO satellites can cover the global regions except for the south and north poles; the signal transmission delay is a constant; frequent inter-satellite switchover is unnecessary; the Doppler shift is small; the technology is matured relatively, and the investment risk is lower, etc. However, a satellite mobile communication system that solely employs GEO satellites has some problems: (1) the GEO satellites have been occupied densely; (2) the south and north poles can't be covered; (3) the signal transmission distance is longer in middle and high latitude regions, since the communication elevation angle thereof is low; (4) high constructions and mountains, etc. between the ground terminal and the satellite hinder the signal transmission, and thus cause shadow regions, so that it is difficult to realize global coverage of satellite mobile communication. For Inclining Geostationary Synchronized Orbit (IGSO) satellites, the sub-satellite track is in a “8” shape with the equator as the symmetry axis, since the orbit inclination angle is greater than 0°; the larger the orbit inclination angle is, the larger the “8”-shaped region is. Therefore, IGSO satellites can effectively overcome one drawback of GEO satellites that the elevation angle in middle and high latitude regions is always low. However, the covering capacity of an IGSO satellite is inferior to that of a GEO satellite. When GEO satellites and IGSO satellites are utilized in combination to complement each other in a network, global covering can be realized.

A multi-beam GEO-IGSO S-LTE mobile communication system employs a large-size antenna array on the satellite to produce multiple beams, which form multiple cells when they reach the ground surface. Similar to T-LTE, GEO-IGSO S-LTE also employs same-frequency networking pattern, in which different main synchronization sequences are configured for adjacent cells. However, GEO-IGSO S-LTE has severe Inter-Beam Interference (IBI), and the overlap range between adjacent cells is large. To implement better cell searching, more main synchronization sequences are required for GEO-IGSO S-LTE.

SUMMARY OF THE INVENTION

Object of the Invention: In view that multi-beam GEO-IGSO S-LTE requires more main synchronization sequences than T-LTE, the object of the present invention is to provide a main synchronization sequence design method for global covering multi-beam satellite LTE.

Technical scheme: To attain the object described above, the present invention employs the following technical scheme: a main synchronization sequence design method for global covering multi-beam satellite LTE, comprising the following steps:

(1) extending a standard Zadoff-Chu sequence into a generalized Zadoff-Chu sequence, so as to obtain more sets of initial candidate main synchronization sequences of candidate sequences;

(2) reducing the sets of candidate main synchronization sequences gradually according to the selection criteria for main synchronization sequences, to obtain a final set of candidate main synchronization sequences;

(3) obtaining a set of main synchronization sequences that has the best trade-off between performance and complexity from the final set of candidate main synchronization sequences.

Furthermore, the step (1) comprises the following steps:

(1) extending a standard Zadoff-Chu sequence in length N_(zc) into a generalized Zadoff-Chu sequence by sub-carrier mapping, to obtain a set of generalized Zadoff-Chu sequences A={x_(u)|u=1, 2, . . . , N_(zc)}, where, x_(u) is the vector form of generalized Zadoff-Chu sequence x_(u)(m), and u represents root index;

(2) removing the middle point from each sequence in the set of generalized Zadoff-Chu sequences A, to obtain an set of initial candidate main synchronization sequences, i.e. B={d_(u)|u=1, 2, . . . , N_(zc)}, where, d_(u) is the vector form of initial candidate main synchronization sequence d_(u)(k), and the length of vector d_(u) is N_(zc)−1.

The generalized Zadoff-Chu sequence x_(u)(m) is obtained from the following expression:

${x_{u}(m)} = {\exp \left( {- \frac{j\; \pi \; {{um}\left( {m + 1} \right)}}{N_{zc}}} \right)}$ m = 0, 1, …  , N_(zc) − 1

where, j is imaginary unit;

${d_{u}(k)} = \left\{ \begin{matrix} {\exp \left( {- \frac{{j\pi}\; {{uk}\left( {k + 1} \right)}}{N_{zc}}} \right)} & {{k = 0},1,\ldots \mspace{14mu},\frac{N_{zc} - 3}{2}} \\ {\exp \left( {- \frac{j\; \pi \; {u\left( {k + 1} \right)}\left( {k + 2} \right)}{N_{zc}}} \right)} & {{k = \frac{N_{zc} - 1}{2}},\frac{N_{zc} + 1}{2},\ldots \mspace{14mu},{N_{zc} - 2}} \end{matrix} \right.$

Furthermore, the selection criteria include: time domain autocorrelation characteristic, peak-to-average ratio of time domain signal, frequency offset susceptibility, time domain cross-correlation characteristic, and computing complexity in the correlation computation of main synchronization sequences on the receiving terminal.

In the step (2), a time domain sequence s_(u)(n) is obtained by carrying out zero padding and N-point inverse discrete Fourier transform for the sequences in the sets of initial candidate main synchronization sequences:

${{s_{u}(n)} = {\frac{1}{\sqrt{N_{zc} - 1}}{\sum\limits_{k = 0}^{N - 1}\; {{{\hat{d}}_{u}(k)}{\exp \left( \frac{j\; 2\; \pi \; {kn}}{N} \right)}}}}},{n = 0},1,\ldots \mspace{14mu},{N - 1}$ ${where},{{{\hat{d}}_{u}(k)} = \left\{ {\begin{matrix} {d_{u}\left( {k + \frac{N_{zc} - 1}{2} - 1} \right)} & {{k = 1},2,\ldots \mspace{14mu},{\frac{N}{2} - 1}} \\ {d_{u}\left( {k - N + \frac{N_{zc} - 1}{2}} \right)} & {{k = {\frac{N}{2} + 1}},{\frac{N}{2} + 2},\ldots \mspace{14mu},{N - 1}} \\ 0 & {{k = 0},\frac{N}{2}} \end{matrix}.} \right.}$

computing the time domain autocorrelation characteristic of the sequences in the sets of initial candidate main synchronization sequences as follows:

$P_{u} = {10\mspace{11mu} {\log_{10}\left( \frac{\max\limits_{\substack{{1 - N} \leq \overset{\sim}{n} \leq {N - 1} \\ \overset{\sim}{n} \neq 0}}\left\{ {R_{u}\left( \overset{\sim}{n} \right)} \right\}}{R_{u}(0)} \right)}\mspace{14mu} ({dB})}$

where,

${{R_{u}\left( \overset{\sim}{n} \right)} = {{\sum\limits_{n = 0}^{N - 1}\; {{s_{u}(n)}{s_{u}^{*}\left( {\left( {n + \overset{\sim}{n}} \right){mod}\mspace{11mu} N} \right)}}}}},$

P_(u) is used to describe the time domain autocorrelation characteristic of a sequence, R_(u)(ñ) is the time domain autocorrelation value of a sequence, * represents conjugation, ñ is the cyclic shift offset of sequence s_(u)(n), and 1−N≦ñ≦N−1;

setting a threshold for the time domain autocorrelation characteristic, and removing the sequences with computed time domain autocorrelation characteristic P_(u) greater than the threshold, to obtain a first set of candidate main synchronization sequences;

computing the peak-to-average ratio of time domain signal of the sequences in the first set of candidate main synchronization sequences:

${{PAPR}\left\{ s_{u} \right\}} = {10\mspace{11mu} {\log_{10}\left( \frac{\max\limits_{0 \leq n \leq {N - 1}}\left\{ {{s_{u}(n)}}^{2} \right\}}{E\left\{ {{s_{u}(n)}}^{2} \right\}} \right)}\mspace{14mu} ({dB})}$

where, PAPR{s_(u)} is used to describe the peak-to-average ratio of time domain signal of a sequence, s_(u) is the vector form of time domain sequence s_(u)(n), E{•} represents statistical average, max{•} represents selecting the maximum value;

setting a threshold for the peak-to-average ratio of time domain signal, and removing the sequences with computed peak-to-average ratio of time domain signal PAPR{s_(u)} greater than the threshold, to obtain a second set of candidate main synchronization sequences;

computing the frequency offset susceptibility of the sequences in the second set of candidate main synchronization sequences:

${{FoS}\left\{ s_{u} \right\}} = {10\mspace{11mu} {\log_{10}\left( \frac{\max\limits_{\substack{{1 - N} \leq d \leq {N - 1} \\ \overset{\sim}{n} \neq 0}}\left\{ {{\overset{\sim}{R}}_{u}\left( \overset{\sim}{n} \right)} \right\}}{{\overset{\sim}{R}}_{u}(0)} \right)}\mspace{14mu} ({dB})}$

where,

${{{\overset{\sim}{R}}_{u}\left( \overset{\sim}{n} \right)} = {{\sum\limits_{n = 0}^{N - 1}\; {{s_{u}(n)}{s_{u}^{*}\left( {\left( {n + \overset{\sim}{n}} \right){mod}\mspace{11mu} N} \right)}{\exp \left( \frac{j\; 2\; \pi \; n\; ɛ}{N} \right)}}}}},$

FoS{s_(u)} is used to describe the frequency offset susceptibility of a sequence, {tilde over (R)}_(u)(ñ) is the time domain autocorrelation value of a sequence after adding a carrier frequency offset, and ε is a normalized carrier frequency offset;

setting a threshold for the frequency offset susceptibility, and removing the sequences with computed frequency offset susceptibility FoS{s_(u)} greater than the threshold, to obtain a third set of candidate main synchronization sequences;

computing the time domain cross-correlation characteristic of the sequences in the third set of candidate main synchronization sequences:

$P_{u_{1}u_{2}} = {10\mspace{11mu} {\log_{10}\left( \frac{\max\limits_{{1 - N} \leq \overset{\sim}{n} \leq {N - 1}}\left\{ {R_{u_{1}u_{2}}\left( \overset{\sim}{n} \right)} \right\}}{R_{u_{1}u_{1}}(0)} \right)}\mspace{14mu} ({dB})}$

where,

${{R_{u_{1}u_{2}}\left( \overset{\sim}{n} \right)} = {{\sum\limits_{n = 0}^{N - 1}\; {{s_{u_{1}}(n)}{s_{u_{2}}^{*}\left( {\left( {n + \overset{\sim}{n}} \right){mod}\mspace{11mu} N} \right)}}}}},$

P_(u) ₁ _(u) ₂ is used to describe the time domain cross-correlation characteristic between sequences, u₁ and u₂ represent different root indexes respectively, and R_(u) ₁ _(u) ₂ (ñ) is the time domain cross-correlation value between sequences;

setting a threshold for the time domain cross-correlation characteristic, and removing the sequence pairs with computed time domain cross-correlation characteristic P_(u) ₁ _(u) ₂ greater than the threshold, to obtain a fourth set of candidate main synchronization sequences;

applying the criteria for computing complexity in correlation computation to the fourth set of candidate main synchronization sequences on the receiving terminal, and utilizing a general sequence property s_(u)(n)=s*_(N) _(zc) _(−u)(n), so that the correlation result of the other sequence can be obtained directly without any additional multiply operation once the receiving terminal obtains the correlation result of one sequence. Thus, a final set of candidate main synchronization sequences is obtained.

Beneficial effects: firstly, global coverage is realized by using GEO satellite and IGSO satellite in combination in a same-frequency network; secondly, with the technical scheme of the present invention, two different sets of main synchronization sequences are provided for GEO satellite and IGSO satellite respectively, so that the interference between cells is reduced compared with cells configured with the same main synchronization sequence, the timing precision is improved, and the cell searching speed is improved; thirdly, with the technical scheme of the present invention, since two different sets of main synchronization sequences are provided for GEO satellite and IGSO satellite, the type of satellite to which a cell belongs can be judged by detecting the main synchronization sequence; fourthly, with the technical scheme of the present invention, multiple main synchronization sequences designed for GEO satellite and IGSO satellite can be two subsets of sequences that are conjugated to each other, and thereby the correlation computation load on the receiving terminal can be reduced.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a flow chart of the main synchronization sequence design method for global covering multi-beam S-LTE disclosed in the present invention;

FIG. 2 illustrates the mapping pattern of a T-LTE main synchronization sequence in the frequency domain;

FIG. 3 illustrates the time domain autocorrelation characteristic of the first set of candidate main synchronization sequences;

FIG. 4 illustrates the peak-to-average ratio of time domain signal of the second set of candidate main synchronization sequences;

FIG. 5 illustrates the frequency offset susceptibility of the third set of candidate main synchronization sequences when the normalized frequency offset is 0.3.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

Hereunder the present invention will be further detailed in embodiments, with reference to the accompanying drawings. It should be appreciated that these embodiments are provided only for describing the present invention, and shall not be deemed as constituting any limitation to the scope of the present invention. Those skilled in the art can make various equivalent modifications to the present invention without departing from the spirit of the present invention, but all these modifications shall be deemed as falling into the protected domain as confined by the attached claims.

The present invention provides a main synchronization sequence design method for global covering multi-beam S-LTE. FIG. 1 is a flow chart of the method. Hereunder the embodiments of the present invention will be further detailed exemplarily in the case of a standard Zadoff-Chu sequence with sequence length N_(zc)=63, with reference to the accompanying drawings.

(1) According to the definition of a standard Zadoff-Chu sequence, N_(zc) and root index u must meet:

gcd(N _(zc) ,u)=1   (expression 1)

where, gcd(x,y) represents the greatest common divisor of x and y. The LTE technical specification specifies that the middle element in the main synchronization sequence corresponding to the DC sub-carrier must be set to zero; thus, a standard Zadoff-Chu sequence can be extended to a generalized Zadoff-Chu sequence, i.e., the value of root index u can be u ∈ {1,2, . . . , N_(zc)}. In that way, more sets of initial candidate main synchronization sequences of candidate sequences can be obtained.

1) Extend a standard Zadoff-Chu sequence into a generalized Zadoff-Chu sequence, to obtain a set of generalized Zadoff-Chu sequences A={x_(u)|u=1,2, . . . , N_(zc)}, where, x_(u) is the vector form of sequence x_(u)(m), and x_(u)(m) is obtained from the following expression:

$\begin{matrix} {{{x_{u}(m)} = {\exp \left( {- \frac{j\; \pi \; {{um}\left( {m + 1} \right)}}{N_{zc}}} \right)}}{{m = 0},1,\ldots \mspace{14mu},{N_{zc} - 1}}} & \left( {{expression}\mspace{14mu} 2} \right) \end{matrix}$

2) Remove the middle point from each sequence in the set of generalized Zadoff-Chu sequences A, to obtain a set of initial candidate main synchronization sequences B={d_(u)|u=1,2, . . . , N_(zc)}, where, d_(u) is the vector form of sequence d_(u)(k), and d_(u)(k) is obtained from the following expression:

$\begin{matrix} {{d_{u}(k)} = \left\{ \begin{matrix} {\exp \left( {- \frac{j\; \pi \; {{uk}\left( {k + 1} \right)}}{N_{zc}}} \right)} & {{k = 0},1,\ldots \mspace{14mu},\frac{N_{zc} - 3}{2}} \\ {\exp \left( {- \frac{j\; \pi \; {u\left( {k + 1} \right)}\left( {k + 2} \right)}{N_{zc}}} \right)} & {{k = \frac{N_{zc} - 1}{2}},\frac{N_{zc} + 1}{2},\ldots \mspace{14mu},{N_{zc} - 2}} \end{matrix} \right.} & \left( {{expression}\mspace{14mu} 3} \right) \end{matrix}$

(2) Reduce the sets of candidate main synchronization sequences gradually according to the selection criteria for main synchronization sequences, i.e., time domain autocorrelation characteristic, peak-to-average ratio of time domain signal, frequency offset susceptibility, time domain cross-correlation characteristic, computing complexity in correlation computation of main synchronization sequences on the receiving terminal, to obtain a final set of candidate main synchronization sequences;

1) Carry out zero padding and IDFT transform with N=64 (i.e., the sampling rate is 0.96 MHz) for the sequences in the set of initial candidate main synchronization sequences B according to the mapping pattern of T-LTE main synchronization sequence in the frequency domain, to obtain a time domain sequence set S={s_(u)|u=1,2, . . . , N_(zc)}, where, s_(u) is the vector form of time domain sequence s_(u)(n), and s_(u)(n) is obtained from the following expression:

$\begin{matrix} {{{{s_{u}(n)} = {\frac{1}{\sqrt{N_{zc} - 1}}{\sum\limits_{k = 0}^{N - 1}\; {{{\hat{d}}_{u}(k)}\exp \left\{ \frac{j\; 2\; \pi \; {kn}}{N} \right\}}}}},{n = 0},1,\ldots \mspace{14mu},{N - 1}}\mspace{20mu} {{where},}} & \left( {{expression}\mspace{14mu} 4} \right) \\ {{{\hat{d}}_{u}(k)} = \left\{ \begin{matrix} {d_{u}\left( {k + \frac{N_{zc} - 1}{2} - 1} \right)} & {{k = 1},2,\ldots \mspace{14mu},{\frac{N}{2} - 1}} \\ {d_{u}\left( {k - N + \frac{N_{zc} - 1}{2}} \right)} & {{k = {\frac{N}{2} + 1}},{\frac{N}{2} + 2},\ldots \mspace{14mu},{N - 1}} \\ 0 & {{k = 0},\frac{N}{2}} \end{matrix} \right.} & \left( {{expression}\mspace{14mu} 5} \right) \end{matrix}$

FIG. 2 illustrates the mapping pattern of a T-LTE main synchronization sequence in the frequency domain.

2) Compute the time domain autocorrelation characteristic P_(u) of the sequences in the set of initial candidate main synchronization sequences, with the following expression:

$\begin{matrix} {{P_{u} = {10\mspace{11mu} {\log_{10}\left( \frac{\max\limits_{\substack{{1 - N} \leq \overset{\sim}{n} \leq {N - 1} \\ \overset{\sim}{n} \neq 0}}\left\{ {R_{u}\left( \overset{\sim}{n} \right)} \right\}}{R_{u}(0)} \right)}}}{{where},}} & \left( {{expression}\mspace{14mu} 6} \right) \\ {{R_{u}\left( \overset{\sim}{n} \right)} = {{\sum\limits_{n = 0}^{N - 1}\; {{s_{u}(n)}{s_{u}^{*}\left( {\left( {n + \overset{\sim}{n}} \right){mod}\mspace{11mu} N} \right)}}}}} & \left( {{expression}\mspace{14mu} 7} \right) \end{matrix}$

Set the threshold for the time domain autocorrelation characteristic as −14.9 dB, and removing the sequences with computed time domain autocorrelation characteristic greater than the threshold, to obtain a first set of candidate main synchronization sequences C={d_(u)|u=1,2, . . . , N_(zc)−1}. FIG. 3 illustrates the time domain autocorrelation characteristic of the first set of initial candidate main synchronization sequences.

3) Compute the peak-to-average ratio of time domain signal for the sequences in the first set of candidate main synchronization sequences obtained in step 2), with the following expression:

$\begin{matrix} {{{PAPR}\left\{ s_{u} \right\}} = {10\mspace{11mu} {\log_{10}\left( \frac{\max\limits_{0 \leq n \leq {N - 1}}\left\{ {{s_{u}(n)}}^{2} \right\}}{E\left\{ {{s_{u}(n)}}^{2} \right\}} \right)}\mspace{14mu} ({dB})}} & \left( {{expression}\mspace{14mu} 8} \right) \end{matrix}$

where, s_(u) is the vector form of time domain sequence s_(u)(n), and E{•} represents statistical average. Set the threshold for peak-to-average ratio of time domain signal as 5.5 dB, and remove the sequences with computed peak-to-average ratio of time domain signal greater than the threshold, to obtain a second set of candidate main synchronization sequences:

$D = \left\{ d_{u} \middle| \begin{matrix} {{u = 1},2,3,4,5,8,10,11,12,13,15,16,19,23,25,} \\ {26,29,30,31,32,33,34,37,38,40,44,47,} \\ {48,50,51,52,53,55,58,59,60,61,62} \end{matrix} \right\}$

FIG. 4 illustrates the peak-to-average ratio of time domain signal of the second set of candidate main synchronization sequences.

4) Compute the frequency offset susceptibility of the sequences in the second set of candidate main synchronization sequences obtained in step 3) with the following expression, where, the normalized carrier frequency offset (abbreviated as “normalized frequency offset”) is set as ε=0.3:

$\begin{matrix} {{{{FoS}\left\{ s_{u} \right\}} = {10\mspace{11mu} {\log_{10}\left( \frac{\max\limits_{\substack{{1 - N} \leq d \leq {N - 1} \\ \overset{\sim}{n} \neq 0}}\left\{ {{\overset{\sim}{R}}_{u}\left( \overset{\sim}{n} \right)} \right\}}{{\overset{\sim}{R}}_{u}(0)} \right)}}}{{where},}} & \left( {{expression}\mspace{14mu} 9} \right) \\ {{{\overset{\sim}{R}}_{u}\left( \overset{\sim}{n} \right)} = {{\sum\limits_{n = 0}^{N - 1}\; {{s_{u}(n)}{s_{u}^{*}\left( {\left( {n + \overset{\sim}{n}} \right){mod}\mspace{11mu} N} \right)}{\exp \left( \frac{j\; 2\; \pi \; n\; ɛ}{N} \right)}}}}} & \left( {{expression}\mspace{14mu} 10} \right) \end{matrix}$

Set the threshold for frequency offset susceptibility as −4.6 dB, and remove the sequences with frequency offset susceptibility greater than the threshold, to obtain a third set of candidate main synchronization sequences:

F={d _(u) |u=19,25,26,29,30,33,34,38,44}

FIG. 5 shows the frequency offset susceptibility of the third set of candidate main synchronization sequences when the normalized frequency offset is set as ε=0.3.

5) Compute the time domain cross-correlation characteristic P_(u) ₁ _(u) ₂ between every two sequences in the third set of candidate main synchronization sequences obtained in step 4), with the following expression:

$\begin{matrix} {{P_{u_{1}u_{2}} = {10\mspace{11mu} {\log_{10}\left( \frac{\max\limits_{{1 - N} \leq \overset{\sim}{n} \leq {N - 1}}\left\{ {R_{u_{1}u_{2}}\left( \overset{\sim}{n} \right)} \right\}}{R_{u_{1}u_{1}}(0)} \right)}\mspace{14mu} ({dB})}}{{where}\text{:}}} & \left( {{expression}\mspace{14mu} 11} \right) \\ {{R_{u_{1}u_{2}}\left( \overset{\sim}{n} \right)} = {{\sum\limits_{n = 0}^{N - 1}\; {{s_{u_{1}}(n)}{s_{u_{2}}^{*}\left( {\left( {n + \overset{\sim}{n}} \right){mod}\mspace{11mu} N} \right)}}}}} & \left( {{expression}\mspace{14mu} 12} \right) \end{matrix}$

Thus, table 1 is obtained.

TABLE 1 Absolute Value of Maximum Cross-Correlation between Sequences (dB) Cross-Correlation( dB) u2 u1 19 25 26 29 30 33 34 38 44 19 0 −6.19 −4.68 −6.30 −6.37 −4.82 −6.24 −6.28 −7.07 25 −6.19 0 −8.06 −7.43 −6.35 −7.02 −4.15 −6.71 −6.28 26 −4.68 −8.06 0 −6.24 −7.43 −4.68 −7.02 −6.24 −4.15 29 −6.30 −7.43 −6.24 0 −8.06 −7.43 −6.35 −4.15 −6.24 30 −6.37 −6.35 −7.43 −8.06 0 −6.24 −7.43 −7.02 −4.82 33 −4.82 −7.02 −4.68 −7.43 −6.24 0 −8.06 −6.35 −6.37 34 −6.24 −4.15 −7.02 −6.35 −7.43 −8.06 0 −7.43 −6.30 38 −6.28 −6.71 −6.24 −4.15 −7.02 −6.35 −7.43 0 −6.19 44 −7.07 −6.28 −4.15 −6.24 −4.82 −6.37 −6.30 −6.19 0

If the threshold for time domain cross-correlation characteristic is set as −6 dB, a fourth set of candidate main synchronization sequences can be obtained by removing the sequence pairs with time domain cross-correlation characteristic greater than the threshold:

G={d _(u) |u=19,25,26,29,30,33,34,38,44}

6) Under the principle of minimum complexity in correlation computation of the main synchronization sequences on the receiving terminal, utilize the following Property:

s _(u)(n)=s* _(N) _(zc) _(−u)(n)   (expression 13)

Set x=x_(R)+jx_(I), y=y_(R)+jy_(I), where, x_(R) and x_(I) are the real part and imaginary part of vector x respectively, and y_(R) and y_(I) are the real part and imaginary part of vector y. Set X_(R)=Circ{x_(R)}, X_(I)=Circ{x_(I)}, where, Circ{x} represent a cyclic matrix in which the first column is x; then, X=Circ{x}=X_(R)+jX_(I), X*=Circ{x*}=X_(R)−jX_(I). The cross-correlation between sequence x and sequence y can be expressed as:

$\begin{matrix} \begin{matrix} {r_{xy} = {X^{H}y}} \\ {= {\left( {X_{R} + {jX}_{I}} \right)^{H}\left( {y_{R} + {jy}_{I}} \right)}} \\ {= {\left( {{X_{R}^{T}y_{R}} + {X_{I}^{T}y_{I}}} \right) + {j\left( {{X_{R}^{T}y_{I}} - {X_{I}^{T}y_{R}}} \right)}}} \\ {= {\left( {a + b} \right) + {j\left( {c - d} \right)}}} \end{matrix} & \left( {{expression}\mspace{14mu} 14} \right) \end{matrix}$

where, a=X_(R) ^(T)y_(R), b=X_(I) ^(T)y_(I), c=X_(R) ^(T)y_(I), d=X_(I) ^(T)y_(R). The cross-correlation between sequence x* and sequence y can be expressed as:

$\begin{matrix} \begin{matrix} {r_{x^{*}y} = {\left( X^{*} \right)^{H}y}} \\ {= {\left( {X_{R} - {jX}_{I}} \right)^{H}\left( {y_{R} + {jy}_{I}} \right)}} \\ {= {\left( {{X_{R}^{T}y_{R}} - {X_{I}^{T}y_{I}}} \right) + {j\left( {{X_{R}^{T}y_{I}} + {X_{I}^{T}y_{R}}} \right)}}} \\ {= {\left( {a - b} \right) + {j\left( {c + d} \right)}}} \end{matrix} & \left( {{expression}\mspace{14mu} 15} \right) \end{matrix}$

Thus, utilizing the conjugation characteristic between the sequences, when the correlation result of one sequence is obtained by the receiving terminal, the correlation result of the other sequence can be obtained directly from the relationship between expression 14 and expression 15 without any additional multiplications. Hence the complexity in correlation computation of the main synchronization sequences on the receiving terminal can be reduced. According to that selection criterion, a final set of candidate main synchronization sequences is obtained:

H={d _(u) |u=(19,44),(25,38),(29,34),(30,33)}

where, a set of 8 sequences is divided into 4 subsets, and the sequences in the sequence pair in each subset are conjugated to each other.

(3) Obtain a set of main synchronization sequences that has the best trade-off between performance and complexity of computation from the final set of candidate main synchronization sequences. To utilize the conjugation characteristic between the sequences to reduce the complexity in correlation computation, 3 sequences can be selected as main synchronization sequences for the GEO satellite, which are selected from any three sequence subsets of the final set of candidate main synchronization sequences H respectively (select one sequence from each of the three sequence subsets), while the IGSO satellite can select the remaining 3 sequences in the 3 subsets. For example, the GEO satellite employs a set of main synchronization sequences H_(GEO)

{d_(u)|u=25,29,30}, while the IGSO satellite employs a set of main synchronization sequences H_(IGSO)

{d_(u)|u=33,34,38}.

It will be understood that changes in the details, materials, steps and arrangements of parts which have been described and illustrated to explain the nature of the invention will occur to and may be made by those skilled in the art upon a reading of this disclosure within the principles and scope of the invention. The foregoing description illustrates the preferred embodiments of the invention; however, concepts, as based upon the description, may be employed in other embodiments without departing from the scope of the invention. The invention is not otherwise limited, except for the recitation of the claims set forth below. 

1. A main synchronization sequence design method for global covering multi-beam satellite LTE, comprising the following steps: (1) extending a standard Zadoff-Chu sequence into a generalized Zadoff-Chu sequence, so as to obtain more sets of initial candidate main synchronization sequences of candidate sequences; (2) reducing the sets of candidate main synchronization sequences gradually according to the selection criteria for main synchronization sequences, to obtain a final set of candidate main synchronization sequences, obtaining a time domain sequence s_(u)(n) by carrying out zero padding and N-point inverse discrete Fourier transform for the sequences in the set of initial candidate main synchronization sequences: ${{s_{u}(n)} = {\frac{1}{\sqrt{N_{zc} - 1}}{\sum\limits_{k = 0}^{N - 1}\; {{{\hat{d}}_{u}(k)}{\exp \left( \frac{j\; 2\; \pi \; k\; n}{N} \right)}}}}},{n = 0},1,\ldots \mspace{14mu},{N - 1}$ ${where},{{{\hat{d}}_{u}(k)} = \left\{ {\begin{matrix} {d_{u}\left( {k + \frac{N_{zc} - 1}{2} - 1} \right)} & {{k = 1},2,\ldots \mspace{14mu},{\frac{N}{2} - 1}} \\ {d_{u}\left( {k - N + \frac{N_{zc} - 1}{2}} \right)} & {{k = {\frac{N}{2} + 1}},{\frac{N}{2} + 2},\ldots \mspace{14mu},{N - 1}} \\ 0 & {{k = 0},\frac{N}{2}} \end{matrix};} \right.}$ (3) obtaining a set of main synchronization sequences that has the best trade-off between performance and complexity from the final set of candidate main synchronization sequences.
 2. The main synchronization sequence design method for global covering multi-beam satellite LTE according to claim 1, wherein, the step (1) comprises the following steps: 1) extending a standard Zadoff-Chu sequence in length N_(zc) into a generalized Zadoff-Chu sequence by sub-carrier mapping, to obtain a set of generalized Zadoff-Chu sequences A={x_(u)|u=1,2, . . . , N_(zc)}, where, x_(u) is the vector form of generalized Zadoff-Chu sequence x_(u)(m), and u represents root index; 2) removing the middle point from each sequence in the set of generalized Zadoff-Chu sequences A, to obtain an set of initial candidate main synchronization sequences B={d_(u)|u=1,2, . . . , N_(zc)}, where, d_(u) is the vector form of initial candidate main synchronization sequence d_(u)(k), and the length of vector d_(u) is N_(zc)−1.
 3. The main synchronization sequence design method for global covering multi-beam satellite LTE according to claim 2, wherein, the generalized Zadoff-Chu sequence x_(u)(m) is obtained from the following expression: ${x_{u}(m)} = {\exp \left( {- \frac{{j\pi}\; {{um}\left( {m + 1} \right)}}{N_{zc}}} \right)}$ m = 0, 1, …  , N_(zc) − 1 where, j is imaginary unit; ${.{d_{u}(k)}} = \left\{ {\begin{matrix} {\exp \left( {- \frac{j\; \pi \; {{uk}\left( {k + 1} \right)}}{N_{zc}}} \right)} & {{k = 0},1,\ldots \mspace{14mu},\frac{N_{zc} - 3}{2}} \\ {\exp \left( {- \frac{j\; \pi \; {u\left( {k + 1} \right)}\left( {k + 2} \right)}{N_{zc}}} \right)} & {{k = \frac{N_{zc} - 1}{2}},\frac{N_{zc} + 1}{2},\ldots \mspace{14mu},{N_{zc} - 2}} \end{matrix}.} \right.$
 4. The main synchronization sequence design method for global covering multi-beam satellite LTE according to claim 3, wherein, the selection criteria include: time domain autocorrelation characteristic, peak-to-average ratio of time domain signal, frequency offset susceptibility, time domain cross-correlation characteristic, and computing complexity in the correlation computation of main synchronization sequences on the receiving terminal.
 5. (canceled)
 6. The main synchronization sequence design method for global covering multi-beam satellite LTE according to claim 4, wherein, the steps further comprising: computing the time domain autocorrelation characteristic of the sequences in the set of initial candidate main synchronization sequences: $P_{u} = {10\mspace{11mu} {\log_{10}\left( \frac{\max\limits_{\substack{{1 - N} \leq \overset{\sim}{n} \leq {N - 1} \\ \overset{\sim}{n} \neq 0}}\left\{ {R_{u}\left( \overset{\sim}{n} \right)} \right\}}{R_{u}(0)} \right)}\mspace{14mu} ({dB})}$ where, ${{R_{u}\left( \overset{\sim}{n} \right)} = {{\sum\limits_{n = 0}^{N - 1}\; {{s_{u}(n)}{s_{u}^{*}\left( {\left( {n + \overset{\sim}{n}} \right){mod}\mspace{11mu} N} \right)}}}}},$ P_(u) is used to describe the time domain autocorrelation characteristic of a sequence, R_(u)(ñ) is the time domain autocorrelation value of a sequence, * represents conjugation, ñ is the cyclic shift offset of sequence s_(u)(n), and 1−N≦ñ≦N−1; setting a threshold for the time domain autocorrelation characteristic, and removing the sequences with computed time domain autocorrelation characteristic P_(u) greater than the threshold, to obtain a first set of candidate main synchronization sequences.
 7. The main synchronization sequence design method for global covering multi-beam satellite LTE according to claim 6, wherein, the steps further comprising: computing the peak-to-average ratio of time domain signal of the sequences in the first set of candidate main synchronization sequences: ${{PAPR}\left\{ s_{u} \right\}} = {10\mspace{11mu} {\log_{10}\left( \frac{\max\limits_{0 \leq n \leq {N - 1}}\left\{ {{s_{u}(n)}}^{2} \right\}}{E\left\{ {{s_{u}(n)}}^{2} \right\}} \right)}({dB})}$ where, PAPR{s_(u)} is used to describe the peak-to-average ratio of time domain signal of a sequence, s_(u) is the vector form of time domain sequence s_(u)(n), E{•} represents statistical average, max{•} represents selecting the maximum value; setting a threshold for the peak-to-average ratio of time domain signal, and removing the sequences with computed peak-to-average ratio of time domain signal PAPR{s_(u)} greater than the threshold, to obtain a second set of candidate main synchronization sequences.
 8. The main synchronization sequence design method for global covering multi-beam satellite LTE according to claim 7, wherein, the steps further comprising: computing the frequency offset susceptibility of the sequences in the second set of candidate main synchronization sequences: ${{FoS}\left\{ s_{u} \right\}} = {10\mspace{11mu} {\log_{10}\left( \frac{\underset{\overset{\sim}{n} \neq 0}{\max\limits_{{1 - N} \leq d \leq {N - 1}}}\left\{ {{\overset{\sim}{R}}_{u}\left( \overset{\sim}{n} \right)} \right\}}{{\overset{\sim}{R}}_{u}(0)} \right)}({dB})}$ where, ${{{\overset{\sim}{R}}_{u}\left( \overset{\sim}{n} \right)} = {{\sum\limits_{n = 0}^{N - 1}{{s_{u}(n)}{s_{u}^{*}\left( {\left( {n + \overset{\sim}{n}} \right){mod}\; N} \right)}{\exp \left( \frac{j\; 2\; \pi \; n\; ɛ}{N} \right)}}}}},$ FoS{s_(u)} is used to describe the frequency offset susceptibility of a sequence, {tilde over (R)}_(u)(ñ) is the time domain autocorrelation value of a sequence after adding a carrier frequency offset, and ε is a normalized carrier frequency offset; setting a threshold for the frequency offset susceptibility, and removing the sequences with computed frequency offset susceptibility FoS{s_(u)} greater than the threshold, to obtain a third set of candidate main synchronization sequences.
 9. The main synchronization sequence design method for global covering multi-beam satellite LTE according to claim 8, wherein, the steps further comprising: computing the time domain cross-correlation characteristic of the sequences in the third set of candidate main synchronization sequences: $P_{u_{1}u_{2}} = {10\mspace{11mu} {\log_{10}\left( \frac{\max\limits_{{1 - N} \leq \overset{\sim}{n} \leq {N - 1}}\left\{ {R_{u_{1}u_{2}}\left( \overset{\sim}{n} \right)} \right\}}{R_{u_{1}u_{2}}(0)} \right)}({dB})}$ where, ${{R_{u_{1}u_{2}}\left( \overset{\sim}{n} \right)} = {{\sum\limits_{n = 0}^{N - 1}{{s_{u_{1}}(n)}{s_{u_{2}}^{*}\left( {\left( {n + \overset{\sim}{n}} \right){mod}\; N} \right)}}}}},$ P_(u) ₁ _(u) ₂ is used to describe the time domain cross-correlation characteristic between sequences, u₁ and u₂ represent different root indexes respectively, and R_(u) ₁ _(u) ₂ (ñ) is the time domain cross-correlation value between sequences; setting a threshold for the time domain cross-correlation characteristic, and removing the sequence pairs with computed time domain cross-correlation characteristic P_(u) ₁ _(u) ₂ greater than the threshold, to obtain a fourth set of candidate main synchronization sequences.
 10. The main synchronization sequence design method for global covering multi-beam satellite LTE according to claim 9, wherein, the steps further comprising: applying the criteria for complexity in correlation computation to the fourth set of main synchronization sequences on the receiving terminal and utilizing a sequence property s_(u)(n)=s*_(N) _(zc) _(−u)(n) to obtain a final set of candidate main synchronization sequences.
 11. A system for designing a main synchronization sequence for a global covering multi-beam satellite LTE having a processor configured to define and generate candidate main synchronization sequences, comprising: a non-transitory computer-readable medium with a sequence of instructions which, when executed by the processor, causes said processor to operate a process to generate and provide candidate main synchronization sequences for a multi-beam satellite LTE, said process comprising: extending a standard Zadoff-Chu sequence into a generalized Zadoff-Chu sequence, so as to obtain more sets of initial candidate main synchronization sequences of candidate sequences; reducing the sets of candidate main synchronization sequences gradually according to the selection criteria for main synchronization sequences, to obtain a final set of candidate main synchronization sequences, obtaining a time domain sequence s_(u)(n) by carrying out zero padding and N-point inverse discrete Fourier transform for the sequences in the set of initial candidate main synchronization sequences: $\underset{\_}{{{s_{u}(n)} = {\frac{1}{\sqrt{N_{zc} - 1}}{\sum\limits_{k = 0}^{N - 1}{{{\hat{d}}_{u}(k)}{\exp \left( \frac{j\; 2\pi \; {kn}}{N} \right)}}}}},{n = 0},1,\ldots \;,{N - 1}}$ ${where},{{{\hat{d}}_{u}(k)} = \left\{ \begin{matrix} {d_{u}\left( {k + \frac{N_{zc} - 1}{2} - 1} \right)} & {{k = 1},2,\ldots \;,{\frac{N}{2} - 1}} \\ {d_{u}\left( {k - N + \frac{N_{zc} - 1}{2}} \right)} & {{k = {\frac{N}{2} + 1}},{\frac{N}{2} + 2},\ldots \;,{{N - 1};{and}}} \\ 0 & {{k = 0},\frac{N}{2}} \end{matrix} \right.}$ obtaining a set of main synchronization sequences that has the best trade-off between performance and complexity from the final set of candidate main synchronization sequences.
 12. The system of claim 11, wherein, the extending step of the process further comprises: extending a standard Zadoff-Chu sequence in length N_(zc) into a generalized Zadoff-Chu sequence by sub-carrier mapping, to obtain a set of generalized Zadoff-Chu sequences A={x_(u)|u=1,2, . . . , N_(zc)}, where, x_(u) is the vector form of generalized Zadoff-Chu sequence x_(u)(m), and u represents root index; and removing the middle point from each sequence in the set of generalized Zadoff-Chu sequences A, to obtain an set of initial candidate main synchronization sequences B={d_(u)|u=1,2, . . . , N_(zc)}, where, d_(u) is the vector form of initial candidate main synchronization sequence d_(u)(k), and the length of vector d_(u) is N_(zc)−1.
 13. The system of claim 12, wherein, the generalized Zadoff-Chu sequence x_(u)(m) is obtained from the following expression: ${{x_{u}(m)} = {{{\exp \left( {- \frac{j\; \pi \; {{um}\left( {m + 1} \right)}}{N_{zc}}} \right)}\mspace{14mu} m} = 0}},1,\ldots \;,{N_{zc} - 1}$ where, j is imaginary unit; ${.{d_{u}(k)}} = \left\{ {\begin{matrix} {\exp \left( {- \frac{j\; \pi \; {{uk}\left( {k + 1} \right)}}{N_{zc}}} \right)} & {{k = 0},1,\ldots \;,\frac{N_{zc} - 3}{2}} \\ {{\exp \left( {- \frac{j\; \pi \; {u\left( {k + 1} \right)}\left( {k + 2} \right)}{N_{zc}}} \right)}\mspace{11mu}} & {{k = \frac{N_{zc} - 1}{2}},\frac{N_{zc} + 1}{2},\ldots \;,{N_{zc} - 2}} \end{matrix}.} \right.$
 14. The system of claim 13, wherein, the selection criteria are selected from one or more of: time domain autocorrelation characteristic, peak-to-average ratio of time domain signal, frequency offset susceptibility, time domain cross-correlation characteristic, and computing complexity in the correlation computation of main synchronization sequences on the receiving terminal.
 15. The system of claim 14 wherein, the process further comprises: computing the time domain autocorrelation characteristic of the sequences in the set of initial candidate main synchronization sequences: $P_{u} = {10\mspace{11mu} {\log_{10}\left( \frac{\underset{\overset{\sim}{n} \neq 0}{\max\limits_{{1 - N} \leq \overset{\sim}{n} \leq {N - 1}}}\left\{ {R_{u}\left( \overset{\sim}{n} \right)} \right\}}{R_{u}(0)} \right)}({dB})}$ where, ${{R_{u}\left( \overset{\sim}{n} \right)} = {{\sum\limits_{n = 0}^{N - 1}{{s_{u}(n)}{s_{u}^{*}\left( {\left( {n + \overset{\sim}{n}} \right){mod}\mspace{11mu} N} \right)}}}}},$ P_(u) being used to describe the time domain autocorrelation characteristic of a sequence, R_(u)(ñ) being the time domain autocorrelation value of a sequence, * representing conjugation, ñ being the cyclic shift offset of sequence s_(u)(n), and 1−N≦ñ≦N−1; and setting a threshold for the time domain autocorrelation characteristic, and removing the sequences with computed time domain autocorrelation characteristic P_(u) greater than the threshold, to obtain a first set of candidate main synchronization sequences.
 16. The system of claim 15, wherein, the process further comprises: computing the peak-to-average ratio of time domain signal of the sequences in the first set of candidate main synchronization sequences: ${{PAPR}\left\{ s_{u} \right\}} = {10\; {\log_{10}\left( \frac{\max\limits_{0 \leq n \leq {N - 1}}\left\{ {{s_{u}(n)}}^{2} \right\}}{E\left\{ {{s_{u}(n)}}^{2} \right\}} \right)}({dB})}$ where, PAPR{s_(u)} is used to describe the peak-to-average ratio of time domain signal of a sequence, s_(u) is the vector form of time domain sequence s_(u)(n), E{•} represents statistical average, max{•} represents selecting the maximum value; and setting a threshold for the peak-to-average ratio of time domain signal, and removing the sequences with computed peak-to-average ratio of time domain signal PAPR{s_(u)} greater than the threshold, to obtain a second set of candidate main synchronization sequences.
 17. The process of claim 16, wherein, the process further comprises: computing the frequency offset susceptibility of the sequences in the second set of candidate main synchronization sequences: ${{FoS}\left\{ s_{u} \right\}} = {10\mspace{11mu} {\log_{10}\left( \frac{\underset{\overset{\sim}{n} \neq 0}{\max\limits_{{1 - N} \leq d \leq {N - 1}}}\left\{ {{\overset{\sim}{R}}_{u}\left( \overset{\sim}{n} \right)} \right\}}{{\overset{\sim}{R}}_{u}(0)} \right)}({dB})}$ where, ${{{\overset{\sim}{R}}_{u}\left( \overset{\sim}{n} \right)} = {{\sum\limits_{n = 0}^{N - 1}{{s_{u}(n)}{s_{u}^{*}\left( {\left( {n + \overset{\sim}{n}} \right){mod}\; N} \right)}{\exp \left( \frac{j\; 2\; \pi \; n\; ɛ}{N} \right)}}}}},$ FoS{s_(u)} is used to describe the frequency offset susceptibility of a sequence, {tilde over (R)}_(u)(ñ) is the time domain autocorrelation value of a sequence after adding a carrier frequency offset, and ε is a normalized carrier frequency offset; and setting a threshold for the frequency offset susceptibility, and removing the sequences with computed frequency offset susceptibility FoS{s_(u)} greater than the threshold, to obtain a third set of candidate main synchronization sequences.
 18. The system of claim 8, wherein, the process further comprises: computing the time domain cross-correlation characteristic of the sequences in the third set of candidate main synchronization sequences: $P_{u_{1}u_{2}} = {10\mspace{11mu} {\log_{10}\left( \frac{\max\limits_{{1 - N} \leq \overset{\sim}{n} \leq {N - 1}}\left\{ {R_{u_{1}u_{2}}\left( \overset{\sim}{n} \right)} \right\}}{R_{u_{1}u_{1}}(0)} \right)}({dB})}$ where, ${{R_{u_{1}u_{2}}\left( \overset{\sim}{n} \right)} = {{\sum\limits_{n = 0}^{N - 1}{{s_{u_{1}}(n)}{s_{u_{2}}^{*}\left( {\left( {n + \overset{\sim}{n}} \right){mod}\; N} \right)}}}}},$ P_(u) ₁ _(u) ₂ is used to describe the time domain cross-correlation characteristic between sequences, u₁ and u₂ represent different root indexes respectively, and R_(u) ₁ _(u) ₂ (ñ) is the time domain cross-correlation value between sequences; and setting a threshold for the time domain cross-correlation characteristic, and removing the sequence pairs with computed time domain cross-correlation characteristic P_(u) ₁ _(u) ₂ greater than the threshold, to obtain a fourth set of candidate main synchronization sequences.
 19. The system of claim 18, wherein, the process further comprises: applying the criteria for complexity in correlation computation to the fourth set of main synchronization sequences on the receiving terminal and utilizing a sequence property s_(u)(n)=s*_(N) _(zc) _(−u)(n) to obtain a final set of candidate main synchronization sequences. 